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TECHNICAL NOTE 2608 

CHARTS AND APPROXIMATE FORMULAS FOR THE ESTIMATION OF 
AEROELASTIC EFFECTS ON THE LOADING OF 
SWEPT AND UNSWEPT WINGS 
By Franklin W. Diederich and Kenneth A. Foss 


SUMMARY 


Charts and approximate formulas are presented for the estimation of 
aeroelastic effects on the spanwise lift distribution, lift-curve slope, 
aerodynamic center, and . damping in roll of swept and .unswept wings at 
subsonic and supersonic speeds. Two types of stiffness distributions are 
considered, one which consists of a variation of the stiffness with the 
fourth power of the chord and one which is based on an idealized constant- 
stress structure. Some design considerations brought out by the results 
of this paper are discussed. 


INTRODUCTION 


A knowledge of the spanwise lift distribution and of some of the 
aerodynamic parameters associated with it is required for the design of . 
a wing structure. Under certain conditions, such as high dynamic pres- 
sures, thin wings, swept wings, or wings designed for low wing loadings, 
the spanwise lift distribution may be affected to a significant extent 
by aeroelastic effects in somewhat the following manner. 

A wing which carries a certain lift necessarily deforms under that 
lift. If the angles of attack along the span are changed as a res ul t of 
this deformation, the lift carried by the wing is changed as well; in 
turn, this change in lift causes a change in the deformation of the wing 
and hence another change in lift, and so on, until an equilibrium condi- 
tion is reached. The changes in the magnitude and the distribution of 
the lift are reflected in changes of the wing lift-curve slope, the wing 
bending and rolling moments, the spanwise center of pressure of the lift, 
and, on a swept wing, the longitudinal center of pressure. Since the 
lift produced by a given change in angle of attack is proportional to the 
dynamic pressure, the various aeroelastic effects tend to increase with 
dynamic pressure. In fact, for certain wings a sufficiently large dynamic 
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pressure may produce a condition of instability in which the change in 
lift caused by deformation is greater than the amount of lift required 
to produce the deformation, so that a given deformation will tend to 
increase until the structure fails. This phenomenon is aeroelastic 
divergence; since it involves only torsional deformations in the case 
of unswept wings, it is often referred to as torsional divergence. 

Several methods are available for calculating these effects (refer- 
ence 1, for instance), but since these effects depend on the structural 
characteristics of the wing, which are not accurately known in advance 
of its design, the relatively large amount of time required for even the 
most efficient of these methods militates against their use in connection 
with preliminary design calculations. A need exists, therefore, for 
means of estimating some of the more important aeroelastic effects on 
the spanwise lift distribution quickly and with an accuracy that is 
sufficient for preliminary design purposes. 

Charts and approximate formulas are presented in this paper for 
estimating the changes in spanwise lift distribution, lift-curve slope, 
wing rolling -moment coefficient, spanwise center of pressure, and aero- 
dynamic center occasioned by aeroelastic action of swept and unswept 
wings at subsonic and supersonic speeds. Also included are summary 
charts which indicate whether the various aeroelastic phenomena con- 
sidered are likely to affect any given design. By means of these charts 
the conventional procedure of designing a wing on the basis of certain 
strength criteria, checking it for aeroelastic phenomena, and then rein- 
forcing it, when necessary, to meet the stiffness requirements imposed 
by these phenomena can often be simplified greatly, inasmuch as the 
effect of some of these phenomena can be estimated in advance of design. 

The use of the charts is described in the section headed "Calculation 
of the Various Aeroelastic. Phenomena, 11 and some considerations involved 
in the selection of the aerodynamic, structural, and geometric parameters 
are discussed in some detail in the section headed "Selection of 
Parameters." These two sections, as well as the sections headed 
"Illustrative Example" and "Preliminary Survey of Aeroelastic Behavior," 
are likely to prove of greatest interest at a first reading of this paper. 
The various parts of the section headed "Discussion" are concerned with 
the limitations of the charts, with the light they shed on such practical 
design problems as the relative significance of strength and stiffness as 
design criteria, with efficient, ways of stiffening a wing that is strong 
but not stiff enough, and with the achievement of aeroisoclinic conditions. 
A brief description of the calculations (based on references 1 and 2) used 
in preparing the charts is contained in the appendixes. 
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SYMBOLS 

aspect ratio (b^/s) 

swept-span aspect, ratio (a/cos^ a) 

cross-sectional area of the (assumed) single torsion cell, 
square inches 

distance from leading edge to section aerodynamic center, 
fraction of chord 

distance from leading edge of mean aerodynamic chord to 
wing aerodynamic center, fraction of mean aerodynamic 
chord 

wing span, inches 

wing span less width of fuselage, inches (b - w) 

wing-root bending-moment coefficient (4M r /qSb) 

lift coefficient of wings alone, exclusive of fuselage 

(Wqs) 

wing lift-curve slope, per radian 

rolling -moment coefficient on both wings alone, exclusive 
of fuselage (Rolling moment/qSb) 

wing-root twisting-moment coefficient ^2T r /qSc r J 

chord (measured perpendicular to elastic axis), inches 

average chord, inches 

section lift -curve slope, per radian 

mean aerodynamic chord (parallel to. plane of symmetry), 
inches 

/ 

Young’s modulus of elasticity, pounds per square inch 
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distance from leading edge to elastic axis, fraction of 
chord 

dimensionless moment arm of the section lift about the 
elastic axis (e - a) 

I 

effective or average dimensionless moment arm 

allowable bending stress, pounds per square inch 
root'- stiffness function 

allowable shear stress, pounds per square inch 
structural weight function 

dimensionless parameters used in approximate formulas for 
angle of attack due to aeroelastic deformation 

dimensionless functions of the distance along the span 
used in approximate formulas for angle of attack due to 
aeroelastic deformation 

modulus of rigidity, pounds per square inch 

wing thickness, inches 

section bending moment of inertia, inches^ 

section moment of inertia in torsion, inches^ 

dimensionless parameters used in approximate formulas for 
dimensionless dynamic pressures at divergence 


dimensionless sweep parameter 


(GJ), 


tan A| 


\e 1 c T (EI) : 

lift of both wings alone, exclusive of fuselage, pounds 
lift per unit distance along span, pounds per inch 


bending moment about an axis perpendicular to elastic axis 
inch -pounds 


free -stream Mach number 
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n 

P 

q 

q* 


q 

s 

s 

s* 

s 

s* 



design load factor 

rolling angular velocity, radians per second 
dynamic pressure, pounds per square foot 

dimensionless dynamic pressure — SI— 

\]M 

dimensionless dynamic pressure I — 2— 

\l44 

total wing area, square inches 

distance along elastic axis measured from wing root, inches 

dimensionless distance along elastic axis (s/s^) 

distance from root to center of pressure of lift along 
elastic axis, inches 

dimensionless distance from root to center of pressure of 
lift (s/s t ) 

accumulated torque about elastic axis, inch-pounds 

thickness of most highly stressed element of skin, inches 

thickness of equivalent skin which includes the material 
in stringers and spar flanges, inches 

distributed torque due to inertia loading, inch-pounds 
per inch 

free -stream velocity, feet per second 
design gross weight of airplane, pounds 
weight of primary structure of both wings, pounds 
weight of both wings exclusive of fuselage, pounds 
width of fuselage, inches 

weight of primary load-carrying structure per unit distance 
along span, pounds per inch 


Ci^ cos A' 
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y* 
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a 


’ll, 2, . . .21 

K 

A 

X 

P 

9 

0 ) 

(i,V,T 

Subscripts: 

cs 

D 


lateral coordinate, inches 

dimensionless lateral coordinate l -2 — \ 

\b/2) 

lateral distance to center of pressure 

angle of attack in a plane parallel to plane of symmetry, 
radians 

angle of local' dihedral, radians, or spanwise slope of 
normal displacement of elastic axis 

density of the material of the primary structure (or an 
equivalent density in the case of sandwich construction) , 
pounds per cubic inch 


lateral distance' measured from wing root, inches 


dimensionless lateral distance 


m: 


V/2 




factors defined in table 1 
ratio of lift-curve slopes 

angle of sweepback at elastic axis 

wing taper ratio (c^/c T ^ 

free-air density, slugs per cubic foot 

angle of structural twist in planes perpendicular to 
elastic axis, radians 

tip stiffness ratio ^(EI^y^E!)-^^ 

dimensionless parameters used in approximate formulas for 
lift, root bending moment, and root twisting moment 

constant -stress 
at divergence 
effective 
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g geometric (due to airplane attitude or built-in twist) 

i inertia 

0 rigid wing (for q = 0) 

r at wing root 

s structural (due to structural or aeroelastic deformation) 

t at wing tip 

Superscripts: 

M due to bending moment 

T due to torque 

f due to root bending 

9 due to root twist 

USE OF THE CHARTS AND A PROXIMATE FORMULAS 
Summary of Method and Scope of Calculations on Which Charts 
and Approximate Formulas Are Based 


Although a detailed understanding of the method and scope of the 
calculations on which the charts of this paper are based is not essential 
to the use of the charts, a brief account of these matters is given, 
primarily to aid in the appreciation of the limitations of the charts. 

The method is described more fully in appendix A. 

Most of the calculations on which the charts are based were made 
by the method of reference 1, which consists in solving the differential 
equations descriptive of an elastically deformed wing under aerodynamic 
loading by numerical methods employing matrix techniques. Treated by 
this method were wings with four taper ratios X (l.O, 0-5, 0.2, and 0), 
two types of stiffness distributions (one proportional to the fourth 
power of the chord and one dictated by constant -stress considerations), 
and four values of a sweep parameter k, at several values of the 

dynamic -pres sure ratio Calculated for each case were the dynamic 

q D • 

pressure at divergence and the changes in spanwise lift distribution, 
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total wing lift, root bending moment, rolling moment, and spanwise center 
of pressure of the lift. For the wings of constant chord and constant 
stiffness, 'calculations were also performed for six’ values of k 'by a 
method which is an extension of that of reference 2 and consists in 
solving the differential equations exactly for these relatively simple 
cases. 

Two important approximations that have been made in all calculations 

are : 


(1) Aerodynamic induction effects at subsonic speeds are taken into 
account by an over-all reduction of the strip-theory loading and, in the 
matrix calculations, by rounding off the strip-theory loading at the tip 
(see references 1 and 2) j at supersonic speeds strip theory is used, with 
a small reduction at the tip for the matrix calculations. 

(2) The. rigid-body rotations imparted to a swept wing by its trian- 
gular root portion are taken into account by a suitable choice of an 
effective root. 

These assumptions were made not so much to simplify the problem as 
to make the results more generally applicable. The most severe limita- 
tion on the use of the charts is probably imposed by the fact that 
calculations have been made for only two types of spanwise distributions 
of bending and torsional stiffnesses: 

(1) Stiffness distributions which vary as the fourth power of the 
chord, such as those of solid wings 

(2) Stiffness distributions associated with structures designed for 
a constant level of combined bending and torsion stress at every point 
on the span, as described in appendix B 

Except for solid wings and those with geometrically similar cross 
sections, for which the stiffnesses vary as the fourth power of the chord 
the stiffness distributions of any given wing depend on the detailed 
design of the wing and cannot be generalized easily. Consequently, the 
constant -stress concept outlined in appendix B has been used to .estimate 
stiffness distributions for some of the calculations; it constitutes an 
effort to relate the stiffness of a wing to its strength on the basis of 
the following assumptions: 

(1) The combined bending and torsional stresses are constant along 
the span. 

(2) The structure is designed for combined bending and torsional 
stresses in such a manner that the sum of the ratio. of the actual to the 
allowable bending stress and the ratio of the actual to the allowable 
torsion stress is equal to unity when the margin of safety is zero. 
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(3) The structure is of the thin-skin, stringer-reinforced shell 
type and its main features do not vary along the span; for instance, the 
number of spars and their chordwise locations are constant along the span. 

(4) At the design condition the spanwise distribution of the applied 
loading is proportional to the chord. 

Although root rotations of swept wings have not been taken into 
account explicitly in the calculations because they vary among different 
designs in a largely unpredictable manner, means for taking these rota- 
tions into account approximately are discussed in a subsequent section. 


Selection of Parameters 

Geometric parameters .- The geometric parameters used in the analysis 
are defined in figure 1. The location of the effective root indicated 
in figure 1 is discussed in the section concerned with the structural 
parameters. 

Aerodynamic parameters .- The aerodynamic parameters which enter the 
analysis are the wing lift-curve slope and the location of the aerodynamic 
center. 'Two lift-curve slopes are used at subsonic speeds: The wing 

lift-curve slope is -used only in conjunction with 'additional lift 

distributions; for all other lift distributions, that is, those due to 
built-in twist, due to roll, or due- to aeroelastic twist, an effective 
lift-curve slope is used. Approximate values of these parameters 

are given for subcritical speeds by the relations 


'hi 


= c- 


A cos A 


-a 


l 


( 1 ) 


■a 


A + 2 — ^ cos A 
2 n 


'V 


= c- 


A cos A 


‘-a 


( 2 ) 


A + b 


1 

t a 

2ji 


cos A 


where c 7 is the lift-curve slope of the section perpendicular to the 
hx 

quarter-chord line at a Mach number equal to M 0 cos A. An approximate 
value is given by 
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c 


2(X 


2* 


/ 


1 - Mq^cos^A 


(3) 


Equation (l) is given in reference 3 and shown to he applicable both to 
incompressible and to subcritical compressible flow. Equation (2) is 
given in references 1 and 2 but without the term c^ 2n in the denomi- 
nator. This term is introduced into equation (2) in order to extend its 
applicability to compressible flows in the same manner as that employed 
for the coefficient of damping in roll in reference 3* 

At supersonic speeds both lift-curve slopes are approximately equal 
to the effective section lift-curve slope 



4 cos A 
/m^cos^A - 1 




The ratio of the lift-curve slopes 




and 


Ct is defined by 
^e 


(5) 


so that for supersonic speeds 


k is equal to 1. 


The local aerodynamic centers are assumed to be at a constant 
fraction of the chord from the leading edge, so that their distances 
from the leading edge (as fractions of the local chords) are all equal 
to the distance of the wing aerodynamic center from the leading edge of 
the mean aerodynamic chord (as a fraction of the mean aerodynamic chord). 
The moment arm e 1 is then given by the relation 

e-^ = e - a (6) 

The lift-curve slopes and the locations of the aerodynamic center 
vary with the free -stream Mach number) hence the appropriate values must 
be used at each flight condition for which aeroelastic calculations are 
made. For airplanes designed to operate at subsonic speeds, only the 
highest Mach number attainable at the highest dynamic pressure is likely 
to be critical from aeroelastic considerations. For airplanes designed 
to operate at supersonic speeds no such general statement can be made; 
however, at a given altitude either the region of Mach numbers near the 
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transition from the subsonic to the transonic regime or the highest 
attainable Mach number is likely to be critical as far as the aeroelastic 
phenomena considered in this paper are concerned. (See figure 4 of 
reference 1 and figure 5 of reference 2, for instance.) 

The airspeeds at which the various aeroelastic phenomena are of 
interest enter the calculations in the form of the corresponding dynamic 
pressures. These dynamic pressures, in turn, are expressed in dimension- 
less form as 


q C I <le e l c r 2s t 2cos A 

^ (gj)_ 


3 * = 


( 7 ) 


or 


q 


a C Iu c rS t 3 sin A 
q e 


U4 (El) r 


( 8 ) 


The parameter q* is useful in the analysis of unswept wings, for which 
torsional deformations are predominant; the parameter "q is useful for 
highly swept wings, for which bending deformations are predominant. In 
general, the parameter q* is used in this paper unless ej_ is zero. 
The ratio of these parameters 


k 


JL - J±- (GJ) r 

q* ' e 1 c r (EI) r 


tan A 


( 9 ) 


is independent of the dynamic pressure and depends only on geometric and 
structural parameters. This' ratio is very useful for analyzing the 
aeroelastic behavior of swept wings. 


Structural parameters . - For the purposes of an aeroelastic analysis, 
the wing structure is characterized by the location of the elastic axis 
and the magnitude and distribution of. the bending and torsional stiffnesses 
(El and GJ), as well as by the magnitude of the rigid-body rotations 
imparted to the wing by its root. 


The elastic axis is usually defined as the locus of points at which 
normal loads. can be applied without causing the- wings to twist. Such a 
locus does not generally exist for practical wings; however, for unswept 
wings without cut-outs an axis can be determined which approximately 
satisfies this condition. Similarly, for swept wings without cut-outs 
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an elastic axis can be defined for the outboard part of the wing if the wing 
is considered to be clamped along some such line as the effective root 
shown in figure 1. In most aeroelastic calculations the locus of shear 
centers for both swept and unswept wings is assumed to be the elastic 
axis. If the structure has large cut-outs which result in sudden changes 
in the stiffnesses and in the shear center along the span, the charts of 
this paper cannot be used except in a qualitative sense. 

The magnitude and the spanwise distributions of the bending and 
torsional stiffnesses enter aeroelastic calculations by means of the 
charts and approximate formulas in different ways. The magnitudes, as 
characterized by the values of the stiffnesses at the effective root, 
have to be known in order to perform any calculations; the distributions 
are implicit in the charts. The root stiffnesses, if not known other- 
wise, can be estimated either from experience with similar designs, from 
the results of the constant-stress concept outlined in appendix B, or 
from a combination of the two. 

The required bending stiffness at the root (EI) r is proportional 

to the design load factor, the wing loading, the wing thickness ratio, the 
fourth power of the root chord, the square of the swept-span aspect ratio 
(a/cos2a), and the ratio of the modulus of elasticity to the allowable 
bending stress, and depends on the taper ratio and on the detailed design 
of the wing (see appendix B) . By means of this relation the bending 
stiffness of one wing can be estimated from that of a reasonably similar 
wing. Or, with the constants of proportionality and F w given in 

table 1 and appendix B, respectively, which take into account some of 
the detailed design parameters as well as the taper ratio, the stiffness 
can be estimated directly. However, in view of the fact that these 
constants have been derived on the basis of a highly idealized structure 
and loading condition they must be used with caution. The ratio of the 
root bending stiffness to the root torsional stiffness can be estimated 
by means of equations given in appendix B or, preferably, from experience 
with structures similar to that under consideration. 

The spanwise stiffness distributions need not be known in detail in 

order to use the charts and approximate formulas. If the wing is solid 

or nearly solid, or if its cross sections are geometrically similar at 
all points, the charts for stiffness distributions proportional to the 
fourth power of the chord are used. If the wing does not have large 

cut-outs and is designed for a constant stress level, the charts for the 

stiffness distributions associated with constant stress are used. The 
use of these charts tends to overestimate aeroelastic effects to some 
extent because, although actual wings are designed for constant stress 
over most of the span, the portions near their tips are designed on the 
basis of other considerations, such as handling loads or minimum standard 
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sheet thicknesses; therefore, wings tend to he stiffer near the tip than 
they would be if designed on the basis of constant stress throughout. 

This difference in stiffness is particularly large if the taper ratio 
is zero* 

If the wing contains large cut-outs or if, for any other reason, 
the wing stiffness distribution is known to be substantially different 
from a constant -stress type, the charts can be used to furnish semi- • 
qualitative results for the various aeroelastic phenomena by using 
fictitious stiffnesses, provided the actual stiffness distribution is 
known at least approximately. The root stiffnesses of these fictitious 
distributiens may be assumed to be the ones that give rise to twist or 
bending angles at the tip which are the same as those of the actual wing 
if the bending moments or tprques vary as the square of the distance from 
the tip. For convenience, ishe spanwise distribution of these fictitious 
stiffnesses may then be assumed to be proportional to the fourth power 
of the chord. On the basis of these assumptions. 


IeiTT = 3X f if* 1 “ s *) 2ds * do) 

e JO 

where the subscript e refers to the fictitious stiffness, and- 

where the integral represents the moment of inertia of the area 

under the function _±_ plotted against s* about the point s* = 1. 

EX 

The fictitious torsional stiffness at the root can he obtained in the 
same manner. The aeroelastic phenomena can then be estimated by use of 
these fictitious root stiffnesses and the charts for stiffness distribu- 
tions proportional to the fourth power of the chord. 

In the derivation of the charts the wing is considered to be clamped 
at the effective root perpendicular to the elastic axis. From the data 
and analyses presented in references 1 , 2, 4, and 5 a satisfactory loca- 
tion of the effective root appears to be at the intersection of the 
elastic axis and the side of the fuselage. 

If the rotations at this effective root are known as a result of 
deflection tests or a detailed analysis such as that of reference 5, the 
root twist due to torque and the root bending due to bending moment may 
be taken into account by moving the effective root inboard by the distance 

AS* = ^( G j) r 


(id 
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or 


r M 

As r = -^(El) r (12) 

where q> is the angle of twist at the root due to a root torque T r , ' 

and where r ^ is the deflection slope at the effective root due to a 
r 

bending moment M r . Since the distances As^ and As r may differ from 
each other, some compromise between the two must be made; for unswept 

wings. the use of As^ appears to be indicated, whereas for highly swept 
wings the use of As^ is more appropriate. 


Preliminary Survey of Aeroelastic Behavior 

The information contained in some of the subsequent sections has 
been summarized in figure 2 for the purpose of ascertaining in advance of 
more detailed estimates, if desired, whether the aeroelastic phenomena 
considered herein are likely to affect the design of the wing structure. 
This preliminary survey is not essential to any of the further calcula- 
tions but may show them to be unnecessary in some cases. 


The charts of figures 2(a) to 2(d) pertain to wings of taper ratios 
0, 0.2, 0.5> and 1.0 and constitute plots of the dynamic -pre s sure param- 
eter q* defined by equation (7) against the sweep parameter k defined 
by equation (9)* These parameters contain the root stiffnesses (GJ) r 

and (El) r ; if, when a preliminary survey of aeroelastic effects is to 
be made, no information whatever concerning the wing stiffness is avail- 
able, the following relations for q* and k may be used: 


q* = 


k = 


C I ^e 1 (l + X) 2 cos A 

i8432 g nW \ „ 

F b S c r ^r 


1 + X G A 


(13) 


( 14 ) 


where F r is a root -stiffness parameter defined in appendix B and 
plotted in figure 3> and where q a and are defined by 
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and 


= ]f_ 

a T l4 T l8 T *9 T ^i5 T 'i7 


(15a) 


% 


^7 

VeVifi 


(15b) 


in terms of some of the factors defined in table 1. 

Figure 2(e) pertains to wings for which the moment arm e 1 is zero 

and, hence, k is infinite; with the degree of approximation involved 
in the use of figures 2(a) to 2(d), figure 2(e) can be used for wings 
with |k|> 2 5. This figure consists in a plot of the dynamic -pres sure 
parameter q, defined by equation (8), against the taper ratio X. If 
no. informat ion_is available concerning the root bending stiffness (El) r 
contained in q, the following relation may be used: 



(l + X)^a^ sin A 


E_ nW ^ ip 

Ft, S c_ n, ’ ^19 F 


r B 


% 


SL - 

3686k 


(16) 


The various lines of the charts of figure 2 designate the conditions 
at which a wing designed on the basis of strength considerations alone 
is likely to encounter divergence or spanwise shifts of the center of 
pressure by various amounts; positive shifts are those toward the tip. 
These spanwise shifts furnish an estimate of the increase in root bending 
moment due to aeroelastic action and an estimate of the shift in wing 
aerodynamic center, since 0 


Aa = sin A -A s 
. c MAC 


(17) 


Inasmuch as the parameters „q* and q contain the dynamic pressure 
the negative values of q* shown in figure 2 may require some explana- 
tion. The four, quadrants of each of the charts of figures 2(a) to 2(d) 
may be characterized for practical purposes as follows: 
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Quadrant 

Sweep 

e l 

Divergence 

Shift in y 

1 

Back 

Positive 

Impossible beyond 
a certain sweep 

Inboard beyond a 
certain sweep 

2 

Forward 

Positive 

Likely 

Outboard 

3 

Back 

Negative 

Impossible 

Inboard 

4 

Forward 

Negative 

Possible beyond a 
certain sweep 

Outboard beyond a 
certain sweep 


For unswept wings k is approximately equal to zero, and the aero- 
elastic phenomena referred to in the charts of figure 2 are similar to 
those of swept wings defined by points in quadrant 2 if q* is positive 
and by points in quadrant '3 if q* is negative. In other words, the 
aeroelastic phenomena of unswept wings are similar to those of swept- 
forward wings if is positive and to those of sweptback wings if 

is negative. The aero dynamic -center shift associated with the shift in 
the lateral center of pressure "y or in the spanwise center of pres- 
sure s is always forward, except for small positive values of k 
(associated with sweep angles smaller than a certain value), in both 
quadrants 1 and k. 


The significance of negative values of q* is that e^ is negative, 
rather than that q is negative. A negative value of e^ may be “obtained 
at supersonic speeds, but under most conditions e^ is likely to be posi- 
tive. Similarly, a negative value of q implies that A is negative 
(that is, that the wing is swept forward). • 


In using figures 2(a) to 2(d), estimates must be made of either the 
root stiffnesses (in conjunction with equations (7) and (9)) or of the 
effectiveness factors q a and (for use in equations (13) and (l4)). 

The factor F r . is obtained from figure 3 for the largest value of e-j_ 
likely to.be encountered at the design load factor and for the given 
taper ratio X. The parameter q* is calculated for the combination of 
dynamic pressure q, lift-curve slope C 


la > 
e 


and moment arm e n which 


is likely to be critical from an aeroelastic point of view. For an 
unswept wing the combination for which the product qC]^ e 1 is a maximum 

is likely to be critical; for a swept wing the combination for which 
qCj^ is a maximum is likely to be critical. The parameter k is then 


calculated for the same value of e 


1 - 
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The values of q*. and k define a point on one of the charts of 
figures 2 (a) to 2 (d) (whichever is closest to the actual taper ratio). 

If the shift in spanwise center of pressure (and any associated shift 
in the aerodynamic center) at that point is small and, in the case of an 
unswept or a swept forward wing, if the absolute value of the ratio of 
the value of q* at that point to the value of q* at divergence for 
the value of k at that point is small, the. static aeroelastic phenomena 
discussed in this paper probably need not be taken into 'account in 
designing the wing structure. On the other hand, if the point on the 
chart indicates the likelihood o‘f significant aeroelastic effects on the 
spanwise center of pressure or the possibility of an approach to the 
divergence condition, further calculations are desirable. The charts 
of this paper may be used for the preliminary calculations} once the 
structure has been designed, more refined methods such as that of refer- 
ence 1 may be used. 

If the moment arm e^ is so small or the angle of sweep so large 
that the parameter k exceeds the range covered by figures 2 (a) to 2 (d), 
the chart of figure 2 (e) may be used for the purpose of a preliminary 
aeroelastic appraisal of the given wing. In this figure only the param- ' 
eter q isjrequired, since k is considered to be infinite. The 
parameter q can be obtained from equation (l 6 ). The analysis then 
proceeds in the same manner as for figures 2 (a) to 2 (d). 


Calculation of the Various Aeroelastic Phenomena 


Dynamic pressure at divergence .- The solutions for the divergence 
speed obtained by the direct method in reference d and those obtained by 
the numerical matrix method given in appendix A can be summarized by 
approximate formulas which give the dimensionless parameters q* D and 

(the values of the parameters defined in equations ( 7 ) and ( 8 ) that 
correspond to the value. of the dynamic pressure q at divergence) in 
terms of their ratio k defined by equation ( 9 ). 

These approximate formulas are • < 




(18) 


and 



( 19 ) 
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When the angle of sweep is zero, equation (l8) reduces to q*^ = K^, 


and when the moment arm e 


^ is zero, as it may he in supersonic flow, 

K, 


— 1 

equation (19) reduces to q^ = - . The constants and 


Kr 


are 


given in table 2 for wings with taper ratios of 0, 0.2, 0.5, and 1.0, 
for both types of stiffness distributions) the parameter q*p for 
unswept wings and the parameter q^ for swept wings with e^ = 0 are 
plotted in figures 4(a) and 4(b), respectively, against the taper ratio 


X. 


With the values of q*p or q^ given by equations (18) and (19) 

and the definitions of these two parameters given by equations (7) 

and (8), the values of q required for divergence may be determined. If 

desired, the corresponding airspeed may be determined from the relation 



The value of q^ is often negative for sweptback wings, and since 

a negative dynamic pressure does not correspond to any real speed, these 
wings cannot diverge. These negative values of q^, nonetheless, are 

useful as reference values in other aeroelastic phenomena. 

The values of the constants and Kg given in reference 2 differ 

somewhat from the corresponding values resulting from the matrix solution 
in appendix A. The matrix results are probably more significant because 
they are based upon more realistic aerodynamic assumptions) the and 

Kg values in reference 2 tend to give conservative results. 

The value of q^ calculated for any given value of q*^ or q^ 
depends on the value of the effective lift-curve slope C T or c 7 

^e CLe 

and, hence, on the Mach number. As suggested in references 1 and 2, the 
value of qp calculated at various Mach numbers may be plotted against 
Mach number. If lines of the actual dynamic pressure at several altitudes 
as a function of Mach number are drawn on the same plot, an intersection 
of the divergence line with one of' the lines of actual dynamic pressure 
designates possible divergence at that value of dynamic pressure, Mach 
number,' and altitude. If this plot is on log-log coordinates, the lines 
of actual dynamic pressure are straight and the ratio of the dynamic 
pressure at divergence to the actual dynamic pressure at a given Mach 
number and altitude can be scaled off directly. (See references 1 and 2.) 
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Spanvise angle -of -attack distribution . - In appendix A, an approxi- 
mate expression is determined for the change in angle of attack due to 
ving flexibility. For the additional -type angle -of -attack distribution 
( a g is constant) the angle of attack due to structural deformation a g 

is given by 


a 

a 


1 

K 


q/qp 

L - -i 



( 20 ) 


The functions and Af 1 , which depend on the spanwise coor- 

dinate s*, and the function F^, which depends on the parameter k, 

are given in figure 5 for swept wings with taper ratios of 0 . 2 , 0 . 5 , 
and 1.0 and with the two different types of stiffness distributions. 

For wings with zero taper ratio the structural deformation cannot be 
obtained from equation (20), as is pointed out in appendix A. However, 
a s 

the ratio k— as a function of the spanwise coordinate s* is shown 

a g 

in figure 6 for the two different stiffness distributions, several values 
of q/q^, and several values of the parameter k. 

The spanwise distribution of a g due to a linear twist (a,g 
which, may be either symmetric or antisymmetric, is approximately 



a 


s 




% 


+ F 2 Af 


( 21 ) 


where the functions f 2 , Af 2 , and F 2 are given in figure 7 for wings 
of taper ratios 0.2, 0.5, and 1.0. The angle -of -attack ratio is shown 
in figure 8 as a function of s*, q/q D , and k for wings of zero taper 
ratio. 


The results of equations (20) and (21) may be superimposed. For 
example, if the spanwise distribution of a g due to rolling is to be 
found, these equations must be added in such proportion that 



y* Ek 
2V 



But 
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so that 


a g = i 1 • V + V s *)il (22) 

where pb/2V is the angle of attack at the tip due to roll. The 


distribution of a s due to roll is then 


a s 

£b 

2V 


q/q D 


-2_ 

q D 


/ 





■-v) 

(fl + P L AfJ + 

2 + F 2 Af 2 ) 


(23) 


Spanwise lift distributions . - If desired, the lift distributions 
can be obtained for the angle-of -attack distributions given in the 
preceding section by one of the commonly used methods of calculating 
spanwise lift distributions, such as that of reference 6. However, the 
following method is simpler. 


■Within the framework of the assumptions made in the analysis the 
lift per inch -of span is proportional to the local angle of attack, so 
that 


7 <Xo 

JL = 1 + K -£ 

V a g 


(24a) 


for geometrical angles of attack which are constant along -the span, and 


Z 


®t 


a, 


s 


-f- = 1 + 

l o V a gt 


(24b) 


a. 


s 


for geometrical angles of attack due to linear twist, where k— and 

g 

<Xg/a g are obtained as indicated in the preceding section. If no better 

approximation is available for the rigid-wing (q = 0) loading 2q, it 
may be estimated from the relation 


V 3 ^ 


( 25 a) 


for geometric angles of attack which are constant along the span, and 
from 


*0 


w cqa 


g 


( 25 b) 


for all other geometric angles of attack. 
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where the coefficients v, n, and t depend on the type of loading. 

The subscript 1 is used for additional -type angle -of -attack distributions 
and the subscript 2 for linear -twist -type angle -of -attack' distributions. 
The coefficients v-^, and are given in figure 9 as functions 

of the parameter k for wings of taper ratios 0.2, 0.5, and 1.0. The 
ratios of the lift, bending -moment, and twisting-moment coefficients to 
their respective rigid-wing values are given in figure 10 as a function 
of q/ q-p for several values of the parameter k when the taper ratio 

is zero. The values of ^2* and^ T 2 are given in figure 11 for 

wings of taper ratios 0.2, 0.5, and 1.0, and ratios of the lift, bending- 
moment, and twisting -moment coefficients are given in figure 12 for wings 
of zero taper ratio. The additional -twist and linear -twist results of 
equations (2 6) to (28) may be superimposed in the same way as those of 
equations (20) and (21). 

The wing rolling -moment coefficient is defined as the rolling 

moment of the loads on both wings about the fuselage center line divided 
by qSb. Therefore, 

2M r cos A + 2T sin A + % K. 

P = ^ W 

*w qSb 


(29) 
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The angle -of -attack distribution due to rolling given in equation (22) 
must be used in finding the values of M r , T r , and L„ in equation (29). 

Spanwise centers of pressure and aerodynamic center . - The spanwise 
location of the center of pressure is given by the distance 


s 



(30) 


or the dimensionless distance 



(31) 


(inasmuch as t|* is equal to s* by virtue of the definitions of those 
dimensionless quantities (see also fig. l), equation (31) can be con- 
sidered to be an expression for q* rather than s*, if desired.) With 
the values of the bending -moment and lift coefficients given in the 
preceding section, the ratio of s to its value for the rigid wing may 
be calculated from either of the equations 


and 


_s 

s’ 


1 - ^L(l - nv) 
q D 


0 1 - -3.(1 - v) 

3D 


As = q D 

8 0 1 - .^1 - v) 


(32) 




where ^ and V-j^ are used for constant geometrical angles of attack 
and' and v 2> for li near ly varying geometrical angles of attack. 

The shift due to aeroelastic action of the longitudinal position of 
the center of pressure, associated with a given shift of the spanwise 
center of pressure. As = s - s Q , is equal to sin A As. The shift in 

aerodynamic center (positive when rearward, or stabilizing) can con- 
sequently be calculated by substituting into equation (17) the values 
of As obtained from equation (32) with values of an ^ v i* 
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Inertia effects .- No charts are presented in this paper for the 
effects of inertia on quasi-static aeroelastic phenomena, that is, aero- 
elastic phenomena associated with flight at constant acceleration; the 
manner in which mass is distributed varies so widely among different 
wings that preparation of a generally applicable set of charts for inertia 
effects appears to be impractical at present. Furthermore, except for 
flying wings, the wing deformations due to inertia loads are small com- 
pared with those due to aerodynamic loads, the two types of loads being 
in about the same ratio as the wing weight to the weight of the entire 
airplane. If desired, however, inertia effects and the aeroelastic 
increment in these effects can be calculated in the manner described in 
the following paragraphs. 

From the known or estimated mass distribution of the wing the 
inertia load Z i per inch of span and the inertia torque t^ per inch 
of span can be calculated for any given normal, pitching, or rolling 
acceleration. Substitution of these loads and torques for the terms 2 
and Ze-^c in equations (A3) or (A 36 ) and equations (A2) or (A35), 
respectively, yields the values of the accumulated bending moment and 
torque due to the distributed inertia loads and torques. In turn, sub- 
stitution of these accumulated bending moments and torques in equa- 
tions (Ak ) , (A 5 ), and (A 6 ), or in equations (A37) and (A 38 ) and the 
matrix equivalent of equation (A 6 ), yields the angle-of-attack distribu- 
tion due to the deformations caused by the inertia effects associated 
with the given acceleration. 

This angle-of-attack distribution can be considered as a geometrical 
angle-of-attack distribution.' For the purpose of calculating the incre- 
ment caused by aeroelastic action, this distribution can be appro xima ted 
by a linear -twist angle-of-attack distribution with a value at the wing 
tip which is such that the moment about the effective wing root of the 
area under the linear -twist distribution equals the moment of the area 
under the calculated angle-of-attack distribution due to inertia effects. 
(The moment, rather than the area, is suggested as a basis of correlation 
because the angles of attack near the wing tip are more important in 
aeroelastic phenomena than those at the wing root.) The justification 
for this rather arbitrary approximation to the angle-of-attack distribu- 
tion is as follows: As previously mentioned, the wing deformations due 

to inertia loads are likely to be small compared with those due to aero- 
dynamic loads; f furthermore , the correction to be applied to these deforma- 
tions as a result of aeroelastic action is usually small compared with 
these deformations and, hence, is very small in comparison with the total 
wing load, so that the correction need not be calculated as accurately 
as .the correction for aeroelastic effects to the rigid-wing lift 
distribution. 

The angle of attack due to structural deformation a associated with 

s 

the lineal-twist distribution can then be obtained from equation ( 21 ) and 



2k 


NACA TN 2608 


figure 7 or, if X = 0, from figure 8. The lift distribution associated 
with the total angle -of -attack distribution due to the deformations caused 
by the inertia effects, including the increment in this angle -of -attack 
distribution produced by aeroelastic action, can then be found from equa- 
tion (2kb) , in which <x g , a g ^, and pertain to the calculated angle- 

of -attack distribution due to the inertia effects (not the linear approxi- 
mation to this distribution) . This lift distribution can be integrated 
to obtain the lift, bending moment, rolling moment, and aerodynamic -center 
position due to inertia effects, as modified by aeroelastic action. 

The lift and rolling moment calculated in this manner may then be 
combined with the lift and rolling moment for steady level or rolling 
flight calculated by the method outlined in the preceding sections. For 
instance, if the contributions of the tail and the fuselage to the air- 
plane lift can be neglected, the wing lift can be written as 


n(W - W w ) 



where 



is the total normal force per unit load factor due to inertia 


effects, including aeroelastic effects; it is equal to -Ww plus the lift 
on both wings due to inertia effects, as modified by aeroelastic action, 
per unit load factor and is almost always negative. In the preceding 
equation Ci^ is a wing, lift-curve slope which includes static aero- 
elastic effects and is equal to C]^ multiplied by the factor on the right 
side of equation (26). Then 


n 


Ct aqS 

1 i ^8 

14k l \ W - W w 

" W - .W w \8n"/s 



144 w - w w 


where 



is a wing lift-curve slope which includes static aeroelastic effects, 
inertia effects, and aeroelastic modification of the inertia effects. 
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Illustrat ive Example 

The parameters of a swept wing, which differs from the wing of the 
illustrative example of reference 1 only in the width of the fuselage to 
which it is attached, are given in table 3- The values of AsT and 
r 

As were calculated from the dimensionless root-rotation constants used 
in the example of reference 1, = 0 and Qp^ = -0.25, by means of 

the relations 

As9 = %p w e 

AsP = Q r M w e 

where w e , as defined in reference 1, is the distance along the span 

between the effective root and the innermost complete section of the 
torsion box perpendicular to the elastic axis. In the wing of the 
illustrative examples of the present paper and of reference 1, w g is 

22. k inches. These relations for As'P and As^ can be obtained from 
equations (ll) and (12) of the present paper in c on j line t ion with the 
definitions of the root-rotation constants given in equations (15a) 
and (l5d) of reference lj in the notation of the present paper the 
definitions are: 

= <Pr T Ar 
w e /(GJ)r 

r m /m 
= r / r 
w e /(El) r 

as the fourth power of the chord in 

the example of the present paper. 

The subsonic and supersonic values of the parameter k were 
calculated from equation (9) • By means of appropriate values of the 
constants and K2 taken from table 3, the values of q*p were 

calculated from equation (l8) and included in table 3* From these values 
of q*pj, the subsonic and supersonic dynamic pressures at divergence were 

found through the use of equation (j) and are given in table 3> These 
values of q D vary as the reciprocal of the effective lift-curve slope, 
the corresponding values of e^ being assumed to remain constant. 



The stiffness is assumed to vary 
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In order to find the angle -of -Attack distribution for additional- 
type loadings from equation (20), the values of F-j_ and of the func- 
tions and Af 1 were taken from figure 5(c). The spanwise change 

in angle of attack is shown in figure 13 for different values of the 
dynamic-pressure ratio. 

The values of t^, and were obtained from figure 9(c) 

and substituted into equations (26), (27), and (28). The wing lift 
coefficient, wing rolling-moment coefficient, and spanwise center-of- 
pressure ratios, as well as the shift in aerodynamic center, were calcu- 
lated by use of these approximate equations in conjunction with equa- 
tions (17) and (29) and are shown in figure 13 as functions of the 
dynamic -pres sure ratio — SL. 


DISCUSSION 

Limitations of the Charts and Approximate Formulas 


The charts and the approximate formulas presented in this paper 
are subject to certain limitations as a result of the approximations 
made in the calculations on which they are based. These limitations 
take the form of restrictions on the plan form, on the speed regime, and 
on the wing structure. The results obtainable by the use of the charts 
are likely to be unsatisfactory for wings of very low aspect ratio or 
very large sweep and relatively unsatisfactory for wings of zero taper 
ratio. 

Wings of low aspect ratio are ruled out on three counts: (l) the 

extent to which aerodynamic forces are overestimated in replacing the 
wing by one with an effective root and tip is larger for wings of low 
aspect ratio than for wings of high aspect ratio, (2) elementary beam 
theory is unsatisfactory for calculating the deformations of wings of 
very low aspect ratio (because the effects of end constraint, shear lag, 
shear deformation, and bending-torsion interaction are more important 
when the aspect ratio is low), and (3) the assumptions made concerning 
the lift distribution of the wing are more nearly true for wings of high 
than for those of low aspect ratio. 

For wings with very large angles of sweep, also, the use of an 
effective root and tip introduces relatively large errors in the aero- 
dynamic forces. Furthermore, the root rotations neglected in the calcula- 
tions (bending rotation due to torsion and twist due to bending) are 
likely to be important for wings with large angles of sweep. 
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The aeroelastic analysis of wings with zero taper ratio entails 
certain mathematical difficulties which do not arise in the case of wings 
with nonzero taper ratio. The stiffness of such wings is zero at the 
tip and very low near the tip, so that the boundary conditions for <p 
and T given by equations (AlOa) to (AlOc) in appendix A are indeter- 
minate. As a result of the relatively large values of the reciprocal of 
the stiffness near the tip, the numerical-integration methods used in the 
matrix calculations are less accurate. These difficulties also occur in 
other methods of solving the aeroelastic equations, such as energy methods. 
Furthermore, the structural behavior near the wing tip is not represented 
adequately by elementary beam theory. Finally, that the aeroelastic 
results calculated for wings of zero taper ratio are not as reliable as 
those for other wings is evidenced also by the fact that they do not 
lend themselves to systematization by means of approximate formulas, as 
do the aeroelastic results calculated for other wings. 


As a result of these considerations one type of plan form of recent 
interest, the delta wing, is seen to be unsuitable for aeroelastic 
analysis by means of these charts because it has a low aspect ratio, a 
large angle of sweep, and zero taper ratio. 

In order to use the charts two aerodynamic parameters must be known 
for any given case, the effective wing lift-curve slope and the section 
aerodynamic center. From an aerodynamic point of view the charts of 
this paper may be used in almost all cases for which these quantities are 
known. The exceptions stem from the fact that the spanwise distribution 
of the lift is assumed to be proportional to the chord, and the distance 
from the section aerodynamic center to the elastic axis (as a fraction 
of the chord) is assumed to be constant along the span. These assumptions 
are not valid for wings with large angles of sweep and wings of low aspect 
ratio, as implied previously. They are also invalid to a greater or 
lesser extent for most wings in the transonic region. Consequently, even 
when the lift -curve slope and the section aerodynamic center are known, 
any results calculated for transonic speeds must be used with caution. 

Another aerodynamic assumption implied in the charts is that no 
concentrated aerodynamic forces, such as those due to a tip tank or 
nacelle, act on the wing. Relatively small nacelies in the inboard half of 
the span can probably be ignored for the purpose of an aeroelastic analysis 
at the preliminary design stage. However, large tip tanks cannot usually 
be ignored even in a preliminary aeroelastic analysis) the aeroelastic 
phenomena may in such cases be greatly underestimated by calculations 
made with the charts of this paper. 


The assumption concerning the applicability of elementary beam 
theory to the calculation of wing deformations due to aeroelastic action 
serves to restrict the wings that can be analyzed by means of the charts 
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to those of moderate or high aspect ratio, as stated previously. Neglect 
of chordwise bending (elastic camber) effects in the calculations on 
which the charts are based serves to impose a lower limit on the thick- 
ness of the wings for which the charts may be used. Whether this limit 
is within the region of practical thicknesses is questionable, however. 
The divergence tests of reference 2, which were performed on flat plates 
of moderately high aspect ratio and with a thickness of 2.5 percent, 
showed no obvious chordwise -bending effects, although the relatively 
small differences between the measured and calculated divergence speeds 
may have been due in part to such effects. 


As mentioned previously, for wings with taper ratios between 0 
and 0.2 the results of aeroelastic calculations are likely to be rela- 
tively unreliable. For taper ratios greater than 0.2, the stiffness of 
actual wings tends to be greater near the tip than that given by the 
constant- stress criterion} consequently, any given aeroelastic effect is 
likely to be somewhat less than that calculated on the basis of a 
constant-stress stiffness distribution, but much larger than that calcu- 
lated on the basis of a c^ distribution. 


If a given structure contains large cut-outs which give rise to 
discontinuities in the stiffness distributions, equation (10) can be 
used to calculate a fictitious root stiffness to be used in conjunction 

with charts for c^-type stiffness distributions, provided the magnitudes 
of the discontinuities are known or estimable. 


Use of the charts of this paper is premised on the assumption that 
the elastic axis is at an approximately constant fraction of the chord. 

If the location of the elastic axis varies somewhat along the span, the 
use of an average value tends to give satisfactory results for the aero- 
elastic phenomena of swept wings} for unswept wings, however, the results 
obtained on the basis of this approximation have to be used with caution. 
If the elastic axis exhibits abrupt ^shifts along the span as a result of 
large cut-outs or for other reasons, the charts should not be used, except 
possibly for moderately or highly swept wings. This restriction is 
mitigated to a certain extent by the fact that ah abrupt shift in the 
locus of shear centers does not necessarily imply an equally large or 
equally abrupt shift in the elastic axis. 


Relation between Strength and Stiffness as Design Criteria 

The strength of a structure is its ability to withstand applied 
loads without failure} the stiffness of a structure is its ability to 
deform relatively little under the applied loads. The two terms are 
related (a fact which forms the basis of the constant -stress type of 
stiffness distributions used in this paper) but are not synonymous. 
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The problem of when to design for strength and when to design for stiff- 
ness and the related problem of how to design a wing for stiffness when 
required to do so have been recognized for a long time. Because of the 
complexity of these problems no generally satisfactory solution exists 
at present, but the charts presented herein shed a certain amount of 
light on the problem insofar as stiffness requirements occasioned by the 
aeroelastic phenomena considered in this paper are concerned. 

The charts of figure 2 indicate the extent to which a wing is likely 
to be affected by aeroelastic phenomena; that is, how far it is from 
divergence and how much its spanwise center of pressure is likely to 
shift as a result of aeroelastic action, provided the wing is designed 
on the basis of strength considerations alone. If the margin against 
divergence is too small, or if the spanwise center of pressure and the 
associated shift in the aerodynamic center are deemed excessive, the 
wing has to be stiffened beyond the amount associated with the required 
strength. The charts of figure 2 therefore serve to delimit the regions 
in which a wing can be designed on the basis of strength considerations 
alone and those in which stiffness considerations predominate, at least 
to the extent of satisfying the stiffness requirements associated with 
the aeroelastic phenomena considered herein. 


The bending moment of inertia required by considerations of strength ' 
alone -for the root section of a wing is directly proportional to the 
design load n(W - W w ), to the spanwise coordinate of the center of pres- 
sure, and to the wing thickness at the ro.ot and is inversely proportional 
to the allowable bending stress Fg. Alternatively, this bending moment 


of inertia may be considered to be proportional to the design wing 
n(W - W ) • 

loading , to the square of the wing area, to the wing thickness 


ratio at the root, and to a function of the taper ratio /which is I 

' (1 + X ) 2 

> this bending moment of inertia is 

inversely proportional to Fg and independent of the aspect ratio. These 
relations for the bending moment of inertia required by considerations of 
strength alone, must be kept in mind in the following discussion of the 
bending moment of inertia required by considerations of stiffness. 


if strip theory is assumed to apply J 


In general, a wing with a high value of q* (see equations (7) 
and (13)) is most likely to be affected by aeroelastic ity (see fig. 2) 
and, for a given value of q*, swept wings are much more likely to be 
affected by aeroelastic ity than unswept ones. (See fig. 2 and equa- 
tions (9) and (14).) Consequently, the following wings are most likely 
to be subject to aeroelastic phenomena, provided they are designed on 
the basis of strength considerations alone: 
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(1) Wings designed for a high flying speed or high dynamic pressure 

(2) Swept wings 

(3) Thin wings 

(k) Wings designed for a low wing loading 

( 5 ) Unswept and moderately swept wings with an elastic axis rela- 
tively far back on the chord or likely to fly in a condition in which 
the section aerodynamic centers are relatively far forward on the chord 

( 6 ) Wings operating at a Mach number at which the lift -curve slope 
is relatively high. 

For given wing loadings and given wing areas, some aeroelastic ' 
phenomena of wings designed on the basis of strength considerations alone 
appear to be substantially unaffected by changes in the taper ratio - for 
instance, the spanwise shift of the center of pressure and the dynamic 
pressure required for divergence, (in the case of the dynamic pressure 
required for divergence, the parameter q*^ (fig. , the root stiffness, 

and the root chord decrease with increasing taper ratio, and the net effect 
of taper is small.) On the other hand, the change in the lift due to 
aeroelastic action is more sensitive to the taper ratio; it is more 
significant for wings with high taper ratio than for wings with low taper 
ratio. 

The effect of aspect ratio on aeroelastic phenomena tends to be 
small for unswept win g s of a given wing area, because these phenomena 
are determined largely by the magnitude of the parameter q*, which is 
independent of the aspect ratio for a given wing area. For the aero- 
elastic phenomena of highly swept wings, however, the parameter q is 
more significant. This parameter is proportional to the swept- span 
aspect ratio for wings of a given area. Consequently, with a given wing 
area, taper ratio, and design wing loading, the aeroelastic effects of 
swept wings tend to be more pronounced for wings with high aspect ratio 
than for those with low aspect ratio. This statement is particularly 
true for the shift of the aerodynamic center, because a given spanwise 
shift of the center of pressure results in a much greater chordwise shift 
in the case of a swept wing of high aspect ratio than in the case of a 
swept wing with low aspect ratio. 


Structural Weight Associated with the Required Stiffness 

When a given wing has been shown to be. subject to undesirably large 
aeroelastic effects (by means of the charts of this paper or by any 
other method), the problem arises how to distribute the additional required 
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stiffness. If, for instance, the pressure on an unswept wing is within 
10 percent of the dynamic pressure required for divergence and a margin 
of 20 percent is desired, an increase of 10 percent in the torsional 
stiffness along the entire span will produce the desired result. The 
question remains, however, whether structural weight can be saved by 
increasing the stiffness more than 10 percent in some places and less in 
others. 


Some insight into this problem may be gained, at least insofar as 
the aeroelastic phenomena considered herein are concerned, from aero- 
elastic and weight calculations that have been made for a family of 
somewhat arbitrarily selected stiffness distributions which- differ from 
the distribution required by the constant -stress criterion. The ratios 
of the local stiffnesses to those associated with constant stress are 
shown at the top of figure 14. The structural-weight factor F w is 
shown for two of these stiffness distributions as a function of the taper 
ratio. The function F w is proportional to the weight W g of the 
primary load-carrying structure and depends on the manner in which the 
wing stiffness and thickness are distributed along the span. (See 
appendix B . ) 


The results of the aerpelastic calculations for wings with taper 
ratio 0.5, constant wing thickness ratio h/c along the span, and these 
two stiffness. .distributions are included in table 2 and figures 5(b), 

7(b), 9(b), and 11(b). The designation "excess strength" in these 
figures refers to the stiffness distributions increased over the constant- 
stress requirement, as shown in figure 14, with a value of oj = 2.0. The 
results of the aeroelastic calculations for the stiffness distributions 
decreased below the constant -stress requirement to a value of <0 = 0.5 
are the same as those for the constant -stress stiffness distributions for 


wings with varying wing thickness ratio; that is, 


(h/c ) t 

(h/ c ) r 


0.5. 


The results of the weight calculations and the aeroelastic calcula- 
tions may be combined in several ways. The dynamic pressure at divergence, 
for instance, can be varied by changing the bending and torsional stiff- 
nesses uniformly along the' span, by leaving the stiffnesses at the root 
unchanged and varying the stiffness distribution in a manner similar to 
that indicated at the top of figure 14, or by a combination of the 
processes. A specified dynamic pressure at divergence can therefore be 
obtained as the result of several combinations of root stiffnesses and 
stiffness distributions. Figure 15(a) consists in. essence of a plot of 
the structural weights associated with various combinations of this 
type against the tip stiffness ratio cs for a specified dynamic pres- 
sure at divergence. This figure indicates that the least weight is 
associated with values of the tip stiffness ratio greater than 1. 
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Similar ly, figures 15(b) and 15(c) consist in essence of plots of the 
structural weights associated with various combinations of root stiff- 
nesses and stiffness distributions required for shifts of ±10 percent 
in the spanwise center of pressure at a specified dynamic pressure. 

^Figures 15(b) and 15(c) also indicate that the structural weight is 
least for values of the stiffness ratio 0) greater than 1. 

The significance of figure 15 is that, if a given wing designed on 
the basis of strength alone needs to be stiffened for aeroelastic reasons, 
most of the stiffening material should be added in the outboard regions, 
provided the weight of the material other than that of the primary load- 
carrying structure is unaffected by the stiffening process. In fact, on 
the basis of aeroelastic considerations alone, weight might be saved in 
some cases by removing material from the root and. adding material at the 
tip) needless to say, however, strength requirements would be violated 
by this procedure . Just where the material should be added in the out- 
board regions cannot be said on the basis of the calculations made for 
figure 15, since these calculations assume any modifications to the 
constant-stress stiffness distributions to be made as indicated at the 
top of figure 14. However, it appears unlikely that great weight savings 
can be had by using modifications which differ substantially from those 
of figure 14. 


The Aeroisoclinic Wing 

The term "aeroisoclinic" refers to wings which deform under an 
aerodynamic load in such a fashion that the angles of attack of all 
sections relative to the free stream remain unchanged. Such a wing has 
the advantage that its aerodynamic loads do not change under aeroelastic 
action either in magnitude or in distribution) its aerodynamic center, 
for instance, is unchanged, and the wing cannot diverge. The achievement 
of such "section aeroisoclinicism" is very difficult and can be realized 
only by separate variation of the bending and torsion stiffnesses) even 
so, the aeroisoclinic condition obtains for only one type of aerodynamic 
loading condition at one Mach number. However, an over-all type of 
aeroisoclinicism in which bending and torsion action tend to cancel for 
the wing as a whole Is relatively easy to achieve. This over-all type has, 
for practical purposes, the same advantages as section aeroisoclinicism, 
in that the aeroelastic phenomena considered in this paper tend to be 
negligibly small for such a wing. 

As may be seen from figure 2, at a small positive value of /the 
parameter k the values of the parameter q* for divergence as well as 
those for given shifts in the spanwise center of pressure tend to infinity. 
This particular value of k represents aeroisoclinic wings in the over- 
all sense) from equations (l8) and (19) it may be seen to be the reciprocal 
of the value of Kq given in table 2. Hence, from the definition of k 
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(equation (9)), 


S t (Gj) r 1 

TeT) 7 tan A = ( 33 ) 


with the implication that the distributions of the stiffness are of 
either the or the constant -s't re ss type and that K 2 pertains to 


Equation (33) 


s t; 


'T> 


either of these types and to the appropriate taper ratio, 
indicates that, for a given plan form with assigned values of 

and A, the disposable parameters for the achievement of aeroisoclinicism 
are the elastic -axis location e^ 

(GJ) r 

and the root -stiffness ratio /■ . — : 

(El) r ' 

the control of the designer. 


which enters into the parameter C p 
the aerodynamic center is not under 


A decrease in the torsional stiffness can sometimes be effected 
without decrease in the bending stiffness or impairment of the strength 
characteristics of the wing, and over-all aeroisoclinicism may be achieved 
in this manner for sweptback wings. Or, if aeroisoclinic, conditions are 
considered -at the outset, a wing can be designed with the elastic-axis 
loeation-relatively far -back . (in the case of a sweptback wing) or forward 
(in front of the aerodynamic center, in the case of a sweptforward wing) 
in order to achieve aeroisoclinicism. However, the fact that only certain 
types of aeroelastic phenomena are considered in this paper must be kept 
in mind. Locating the elastic axis far back or decreasing the torsio nal 
stiffness, for instance, may lead to flutter difficulties, the solution 
of which may require excessive mass balancing of the wing as a whole. 


Relation of Charts to Design Procedure 

The first step in the design of a wing structure, once the wing 
geometry and the over-all airplane characteristics have been decided 
upon, usually consists of a rough apportioning of structural material 
along the span in a manner intended to satisfy strength requirements 
approximately. At a later stage in the design procedure the structure 
is checked for aeroelastic effects and modified, if necessary. The 
modifications are then checked again, and so on, until both stiffness 
and strength requirements are met with what is believed to be a near- 
optimum structure from weight considerations. The charts of this paper 
may be used to facilitate the procedure at several stages. 

At the very outset, the preliminary- survey charts can be used to 
establish some over-all aeroelastic characteristics of the wing structure 
that would be obtained by designing the wing for strength alone. If 
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these characteristics are satisfactory, the design of the wing structure 
can proceed on the basis of strength requirements alone. The final 
design can then be checked for the aeroelastic effects considered in this 
paper by means of the charts contained herein, and for other aeroelastic 
effects, such as flutter and loss of lateral control, by equally approxi- 
mate methods. However, if the preliminary survey indicates that a wing 
designed on the basis of strength alone would be unsatisfactory from 
consideration of aeroelasticity, sufficient additional stiffness may be 
incorporated in the preliminary design stage, provided the taper ratio 
does not differ greatly from 0.5 and the wing thickness ratio is constant 
along the span. For instance, the preliminary- survey charts may indicate 
a shift in the spanwise center of pressure which gives rise to a shift 
of 4 percent in the aerodynamic center, whereas the desired maximum shift; 
is 2 percent, so that the spanwise shift must be reduced to 50 percent of 
that indicated on the preliminary- survey chart. The shifts in the span- 
wise center of pressure for a wing with increased stiffness at the tip 
(the "excess strength" case, for which <jd = 2.0) and for a wing with 
decreased stiffness at the tip (the wing with to = 0.5, for which the 

(h/c) t \ 

results of the case of =0.5 may be used J can then be obtained 

from figure 9(b) and equation (32), in conjunction with the value of the 
dynamic pressure at divergence estimated from equation (l8) or (19) • The 
fact that the wings with o> = 2.0 and o> = 0.5 have different dynamic 
pressures at divergence than does the constant -stress wing must be kept 
in mind. 


From the stiifts of the spanwise center of pressure calculated in 
this manner the value of o> for the desired spanwise shift can be 
obtained by interpolation and, hence, the approximate magnification 
factors to be applied to the stiffness distribution for constant stress 
can be obtained from the chart at the top of figure l4. Estimates for 
the other aeroelastic characteristics considered in this paper can then 
be obtained for the wing with this modified stiffness distribution by 
interpolating between the results given for these characteristics for 
wings with cu = 0.5, 1.0, and 2.0; that is, for the cases referred to, 
respectively, as 


(h/c) t 

(h/c) r 


0.5 


(h/c) t 

r , T ' r = 1.0 

(h/c) r 

and 

(h/c) t _ 
(h/c) r " 


1.0 (excess strength) 
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in table 2 and figures 5(h), 7(h), 9(h), and 11(b). Once the structure 
of such a wing has been designed, the various aeroelastic effects con- 
sidered herein should be checked by a more accurate method, such as that 
of reference 1, and the loss of lateral control and the flutter character- 
istics should be calculated. 


CONCLUDING REMARKS 


Charts have been presented for the estimation of aeroelastic effects 
on the spanwise lift distribution, lift-curve slope, aerodynamic center, 
and damping in roll of swept and unswept wings at subsonic and supersonic 
speeds. Two types of stiffness distributions have been considered, one 
which consists of a variation of the stiffness with the fourth power of 
the chord and is appropriate for solid wings, and one which is based on 
an idealized constant -stress structure and is believed to be more nearly 
representative of actual structures. 

" The limitations of these charts are that they do not apply to wings 
with very low aspect ratio or very large angles of. sweep,, nor to wings 
with large sources of concentrated aerodynamic forces. T he charts are 
likely to be less reliable for wings with zero taper ratio than for wings 
with other taper ratios, and less reliable when the component of the 
Mach number perpendicular to the leading edge is transonic than when this 
component is either subsonic or supersonic. Wings with large discon- 
tinuities in the spanwise distribution of the - bending or torsional stiff- 
nesses cannot be analyzed directly by use of the charts, but a means of 
making approximate calculations for such wings has been presented. No 
charts have been presented for inertia effects but a method of estimating 
these effects has been outlined. 

In addition to facilitating the calculation of various static aero- 
elastic phenomena, the charts serve to simplify design procedure in many 
instances, because they can be used at the preliminary design stage to 
estimate the amount of additional material required to stiffen a wing 
which is strong enough and because they indicate that the best way of 
distributing this additional material is to locate most of it near the 
wing tip. 

Also, the charts facilitate the achievement of aeroisoclinic condi- 
tions, inasmuch as they serve to define a simple relation between the 
elastic -axis location and the wing stiffness ratio which is required to 
obtain this condition for a given plan form. Finally, the charts indicate 
that a wing which is strong enough is most likely to be affected by aero- 
elastic phenomena if it is to operate at high dynamic pressures, if it 
is thin, if it has a large angle of sweep, if it is designed for a low 
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wing loading, if it has an elastic -axis location relatively far hack on 
the chord, and if it is to operate at transonic or high supersonic Mach 
numbers. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., September 13, 1951 
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APPENDIX A 

METHODS OF CAPSULATIONS ON WHICH CHARTS ARE B ASED 
The Aeroelastic Equations 


The methods of calculating aeroelastic phenomena used in preparing 
the charts of this paper are based on the following assumptions: 

(1) Aerodynamic induction is taken into account by applying an 
over-all correction to strip theory and, when matrix integrations are 
used, by rounding off the resulting load distribution at the tip. 

(2) Aerodynamic and elastic forces are based upon the assumption of 
small deflections. 

(3) The wing is clamped at the root perpendicular to a straight 
elastic axis (see fig. 1), and all deformations are considered to be 
given by the elementary theories of bending and torsion about the elastic 
axis . 

In keeping with assumptions (l) and (2), the force per unit width oh 
a wing section perpendicular to the elastic axis is 


1 = + c i tLe a ' s ) 


(Al) 


where a 8 and a g are, respectively-, the angle of attack due to 
structural deformations and the rigid-wing angle of attack, in planes 


parallel to the plane of symmetry. (The geometrical angle of attack is 


considered to be constant along the span in equation (Al)$ 
of linear twist the coefficient C T is used instead of 


in the case 

c v) The 


torque of this force about the elastic axis is Ze-.c for uncambered 
sections . * - 1 


The integral equations for the accumulated torque and the bending 
moment are ~ 


T = 



Ze^c ds 


(A2) 
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(A3) 


and, insofar as assumption (3) holds, the angles of structural twist and 
bending referred to axes parallel or perpendicular to the elastic axis 
are 


<P 


r 



The angle of attack due to structural deformations is related to cp 
and f by the equation 


a s = cp cos A - T sin A 


(A6) 


Combining equations (Al) to (A6) gives two simultaneous differential 
equations: 



2 

" qe l c r 

iior 1 c i u a g + c iu (cp C0B A ■ - r sin A) 


(A7). 


fl( EI IS’ fk 1 ' * ‘V’ C0S A - r Sln A U <A8> 

as 

These equations are subject to the following boundary conditions: 
Zero twist and bending at the root, 

9(0) = 0 (A9a) 

(A9b) 


r(o) = 0 
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Zero torque , moment, and shear at the tip, 



(AlOa) 


(AlOb) 
( AlOc ) 


In the following sections, equations (A7) and (A8) are solved 
explicitly for an untapered wing with constant stiffness along its span 
and by matrix integration for a wing with any arbitrary stiffness and 
chord variation. 


Solutions for Uniform Wings 


Arbitrary geometric angle of attack .- If the torsional stiffness, 
the bending -stiffness, and the chord of the wing have constant 
values (GJ) r , (El) r , and c r , respectively, along the wing span, 

equations (A7) and (A8) become 


cp" cos A = 


Iff 

f sin A = 



+ (cp cos A - f sin A) 


V) 


+ (cp cos A - f sin A) 


1 


(All) 


( A12) 


where the differentiation denoted by the primes is with respect to 
t s 1 - and the dimensionless parameters q* and q are defined by 


C]^ qe^c r 2 s^.2 C os A 


3* = 


iMKgj), 


(A13) 


Cr tt e qcrs t 3sin A 
q lU(El) r 


(■AHO 
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Differentiating equation (All) once with respect to I and combining 
it with equation (A12) yields the single differential equation 


a'” + q*a, 


- = 


K It 


(A 15 ) 


(The factor ^ is used with a for the sake of consistency, despite 

the fact that a geometrical angle of Attack which is constant over the 
span does not have a spanwise derivative.) Equation (A15) is subject to 
the ‘ following boundary conditions: 


From equations (A9a) and (A9b) 

a s (l) = 0 

From equations' (AlOa) and (AlOb) 

a B '(0) = 0 

From equations (AlOc) and (All) 

a s "(0) = -q* 


. ( 0 ) 

— + a s (o) 


(Al6) 


(A17) 


(A18) 


where functional notation is used, . so that, for instance, a s (D means 
the value of a g at I = 1. 

The solution of equation (A15) can be effected very readily by means of 
of Laplace transforms. The complete solution of this equation is 

a s ( ^) =^T) f 3 (l) - H (0 (A19) 

where the integral H(|) is defined as 



q*f 4 (l - l x ) - qf 5 (l 



a gUl) 


d!i 


H(|) = 


K 


(A20) 
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The functions f^, and are defined by 


f 3 (l) =C ie ' 2Pl + e Pl 


C 2 cos ?£ + 22 sin 7^ 


f^Cl) -C^e -2 ^ + e^^C^ cos 7 I + — sin 7i) (A21b) 

f 5 (0 = Cr J e 20 ^ + e0 ^( c 8 cos 7^ + — sin 7 ^ (A21c ) 

-2(3 and M i7 are the roots of r^ + q*r - <f = 0 and 

c = ^P 2 

9P 2 + 7 2 


Co = 1 - C 


c = 3P 3 - P7 2 
3 9P 2 + 7 2 


C, = — -23 
9P 2 + 7^ 


C 5 ~ -C 4 


33 2 + 7 2 
9P 2 + 7^ 


C7 = ^T7 


c 8 — c 7 


Cg = 

° qr2 . ~2 
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The condition for divergence is that a s he finite when a is 

o 

zero along the entire span. As can be seen from equations (A19) and. (A20), 
divergence can occur only when 


f 3 (lj = 0 (A22) 

Thus, for a particular value of the parameter k = ^ the value of q* 
(or q") at divergence is the one which satisfies equation' (A22) . 

Constant geometric angle of attack .- For the additional- type angle 
of attack, a g (|) = Constant: 

H(|) = [l - f 3 (l)]^ (A23) 


and 

* a 3 ( iX + a g _ ££ I) 

a g = f 3 (l) 


(A24) 


For lift distributions ba.sed on assumption (l) given at the beginning 
of this appendix, the lift per unit width of span may then be written as 

_L = V £ 

The wing lift coefficient, the wing-root bending-moment coefficient, 
and the wing-root twisting -moment coefficient are given in general by 


C 


^ <J S 


= c 


s t c r 
L a S 72 “ 



c 




+ a g( u] d! 


(A25) 


4^ 

n _ qSb 


= C 


la 


+a6( ^ as ** 


(A26) 
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C m = 


2T 


T qSc- 


= c io;s7r e i J 0 + “g (5 0 aE 

Then, for the uniform wing. 


(A27) 


\ r 1 Ka s (^ + a , 


JO 


a n 




f^(i) 

fjnr 


(A28) 


C B _ f 1 P 5 - *a B (!) + otg 
C~ = 2 — ~ d S d t 


B 0 JO JO 


g . 


f 5(l) 


" 2 fTHT 


(A29) 


and 


C T 0 = V 


Linearly varying geometric angle of attack .- For the linear-twist- 
type angle of attack, a g (l) = (1 - IJa^, the factor k is 1, and 

H(5) =a g (|) - p 3 (^) - f] f (l)|a gt (A30) 


so that 


l _ a s ( £) + tig( £) 

ST % 

fkd) 

~ f 3 (i) f 3^ " 


(A31) 


The ratios of the wing lift, wing-root twisting moment, and wing- 
root bending-moment coefficients to their rigid-wing values are then, 
on the basis of assumption (l). 
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% _ 2 x ^gU) 

C w 0 Jo a gt 


! f*(l) - f 5 (l) 


c B_ = 3 P P a sU) + a gU) 

c b q Jo Jo a g t 


= 3 ]fW f 5 (l) ' lE3 (l) + q * f 5 (l) " ^ 


lliidl dg 


Ct _ c Iy 
% C Ivo 


as in the preceding section. 


Solution for Nonuniform Wings 

Equation (Al) may be written in the matrix notation of reference 1 
as 

« - ilir ^L'lf s ♦ ' <* 

and equations (A2) and (A3) as 

{T} = Bt[l ’]{le x c} 

= TO C I^ e e l c r 2a t[ r ] jf4 + *{ a g}' (A: 

\ 

{ M } = «t s [ji3 0) 


" lW C I tle c r B t 2 [ IZ 3 ■ {%} + l{ a g) 
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where the matrix Cl '3 performs an integration of the running torque 
Ze^c from the tip inboard, and the matrix [jl3 performs a double 
integration of the running load from the tip inboard. These matrices 
are derived and given in reference 1. They are based upon Simpson's 
with a modification at the tip, where the load distribution is 
assumed to go to zero with an infinite slope at the tip. 

Equations (A4) and (A5), written in matrix notation, are 


= TEI57&I 


(GJ) r 
_ GJ 

(El) r ~ 

El 


(t) 

(A37) 

{M} 

(A38) 


where the matrix ClJ" serves to integrate the accumulated torque or 
bending moment outboard from the wing root. This integrating matrix is 
based upon Simpson’s rule without the tip modification and is given in 
reference 1. 


The substitution of equations (A35), (A 3 6 ), (A37), and (A 38 ) in the 
matrix equivalent of equation (a 6 ) yields 

N = ^[a](h + - (A39) 

where the aeroelastic matrix [aJ is defined by 


DO * H 


The parameters a* and q are defined by equations (A13) and (Al4), 
respectively, and ’ 



k = -2_ 

q* 

When a g is zero along the entire span, equation (A39) becomes 

{%} = q*CA]{a s j • (A4l) 
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Consequently, for a particular value of k, the value of q* at 
divergence can be found by the iteration of the aeroelastic matrix ^Aj* 


Equation (A39) may be rearranged as follows: 

[Lil - q*C A l]{ Ka s + a g} = { a g} 


(A42) 


The set of linear simultaneous equations represented by equation (a 42) 
may then be solved for the total angle of attack Ka g + a g in terms of 
the values of a g along the span. 

The integrations in equations (A25), (A2 6 ), and (A27) may be per- 
formed with the first rows of the jjE'J and [XI 0 matrices. Thus 






r o 


L 1 ' J lLc r ]f a s + a g} 

WiLilN 


(A43) 


_B_ 

CB 0 


Mil??] 

{ Ka s + a g} 

IP'Jl 

C 

c r 

N 


(a44) 


and 


■'Tq 


u'jib 

k)% 

[ Ktx s + a g} 

Lrj 1 

OJ 

1 


(a45) 


The aeroelastic characteristics of uniform wings were calculated by 
both the direct method of the preceding section and the matrix method 
given in this section. The values of the divergence parameter q* D , 

calculated by the direct method, were found to be about 5 percent greater 
than the corresponding values calculated by the matrix method. This 
discrepancy can be shown to be almost entirely due to the rounding off 
of the loading of the wing tip in the matrix method. The differences 

a s ^Lw Cg Cij 

between corresponding values of — , ^ , p — , and — are 

a g c l wq c Bo Lt o 

negligible. 



Combination of Results 


The forms of the approximate formulas used in combining the results 
of the many computations indicated in the analysis were obtained by 
considering a highly idealized semirigid wing; that is, a wing which is 
rigid along its. entire span but can bend and twist at the wing root 
subject to the restraint of a bending and a torsion spring. 


If it is assumed that the two spring constants correspond to 

, (El) r 

and — , the value of q* at divergence is given by the simple 
formula 


(GJ) r 

~ E t~ 



(Ak 6) 


where the factors K-j. and Kg depend on the taper ratio and the spanwise 

variation of the stiffness. As shown in reference 2, this formula serves 
as a good approximation to the calculated values of q* D . 


For the semirigid wing, the ratio of cc s to a.g is found to be 
proportional to j n order to adapt this expression to the 

1 - JL 

flexible wings considered in the present analysis, the following approximate 
expression was found to provide satisfactory correlation: 


a s i q/q D , 

5— = K 7TlL ( + F (A47) 

where f and Af are functions of the spanwise coordinate s* and the 
wing-chord and stiffness variations; F is a function of the parameter k 
and the wing-chord and stiffness variations. The functions f, Af 
and F also depend on the type of spanwise variation of the geometrical 
an €>l e of attack, the subscripts 1 and 2 being used to differentiate 
between the two types of interest. The accuracy of equation (AV7) is 
illustrated in figures l6 and 17. 

If equation (A25) is used for the wing lift coefficient ('with S 
replaced by s*) and equation (AU7) for the angle -of -attack distribution, 
the wing lift coefficient may be expressed as 
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C ^v = v 


2p 

— — C T + C T 

- - 2 - ^0 


(A 48) 


or 


V 


1--JL 

2 D 


• (A49) 


where the parameter 




t — (f + F Af)ds*. 
c r 


So- Qf. 

t — ds* 


JO 


c r % 


(A50) 


so that • v and 


'^Vr 


are functions of k and of the wing-chord and 


stiffness variations and depend on the type of geometrical angle-of- 
attack distribution as well. 


As indicated by equation (a 47), within the approximation inherent 
in that equation, the shape of the spanwise distribution of a s does 
not vary with dynamic pressure. Therefore, to a good approximation, the 
lateral center of pressure of the lift due to a s (as well as that due 
to a,g) does not change its position along the elastic axis when the 
dynamic pressure changes. The following approximate formula for the 
wing-root bending -moment coefficient may therefore be deduced from 
equation (A48): 



s s* v 


_q_ 

2d 


1 - 77 ^ 


'^W , r 


+ s *C T 

g L 'Wf 


(A51) 
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where s g * and s g * are the dimensionless moment arms about the effec- 
tive wing root of the lifts due to a s and a and are defined by 






2 


and 




s * = 




8 i ^ H 
2 t 


Then, 


Cg 1 - ^ (1 ‘ ^ 


-'B, 


0 


1 - -a. 

q D 


(A52) 


where u is defined by 


H-ii* 

§ g* 


(A53) 


so that is a function of the parameter k, of the taper ratio, 

and of the stiffness distributions; it also depends on the type of 
geometrical angle -of -attack distribution. 

An approximate formula for the wing-root twisting-moment coefficient 
may be deduced from equation (A48) as follows: 


C rp — 




_q_ 

q D 


Cj + 6-1 Cy 

- -2- • ^0 lg N> 


(A54) 
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where e-. and e-i are the effective dimensionless moment arms about 
1 s x g 

the elastic axis of the lifts due to a s and <x g and are defined by 


and 






r 


Then 


where 



1 - ^-(1 - vt ) 
q D 

1 - -a. 

q D 


(A55) 



(A56) 


so that is a function of k, the taper ratio, and the stiffness 

C T 0 

distributions and also depends on the type of geometrical angle -of-attack 
distribution. 


The values of v, n, and r are given for the two types of 
geometrical angle -of-attack distributions in figures 9 and 11. 

Figure 18 shows the approximate formulas (Ak9), (A52), and (A55) to 
be in good agreement with more accurately computed values. 

The foregoing approximate formulas for the structural angle of attack 
and for the lift, bending-moment,’ and twisting-moment coefficients are not 
applicable to wings with zero taper ratio. An attempt was made to combine 
and systematize the results calculated for such wings in the manner 
employed for wings with other taper ratios, but the approximate formulas 
obtained in this way were found to yield unreliable results. Consequently 
they are not presented in this paper) instead, the results calculated for 
the wings with zero taper ratio are presented directly in figures 6, 8, 

10, and 12 . 
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APPENDIX B 

STIFFNESS DISTRIBUTION OF CONSTANT-STRESS WINGS 
Outline of Constant -Stress Concept 


In order to calculate aeroelastic effects, the bending and torsional 
stiffnesses of the wing structure, El and GJ, have to be known. 

These stiffnesses enter the calculations in two ways. The root stiff- 
nesses, as indices of the over-all bending and torsional stiffnesses, 
constitute primary parameters which are required for use of the charts 
of this paper, but were not required in the preparation of the charts. 

On the other hand, the stiffness distributions,' that is, the ratios of 
the local stiffnesses along the span to the root stiffnesses, are 
secondary parameters which are not required for use of the charts but 
did have to be assumed in order to prepare them. 

In calculations preliminary to the actual design of the structure, 
the bending and torsional stiffnesses of the structure are not known; 
they must be estimated on the basis of either past experience or considera- 
tions of an idealized structure. For the purpose of estimating stiffness 
distributions, past experience with similar structures is likely to be a 
useful guide in any specific case but does not lend itself to generaliza- 
tion and hence to the preparation of generally applicable charts. The 
stiffness distributions (other than those which vary as the fourth power 
of the chord) used to prepare the charts of this paper have been obtained 
from considerations of an idealized structure, as outlined in this appendix. 

Basically, the method of this appendix consists in an effort to relate 
the stiffness of a wing to its strength and to estimate that strength on 
the basis of certain assumptions. The fundamental assumptions are that 
the bending and torsional stresses are constant along the span and that 
the applied loading is proportional to the local chord. The other assump- 
-tions concern the bending and torsional stresses caused by this load and 
their relation to their allowable values. In estimating these stresses 
the structure. is assumed to be of the thin-skin, stringer-reinforced shell 
type. Certain effectiveness factors are used; for instance, the ratio of 
the allowable torsional stress to the allowable bending stress, or the 
ratio of the cross-sectional area of the effective torsion cell to the 
product of the chord and the wing thickness. The root stiffnesses 
estimated by the method of this appendix depend directly on the values of 
these ratios. The stiffness distributions, on the other hand, are largely 
independent of these ratios but imply the assumption that the ratios are 
approximately constant along the span. Consequently, the constant -stress 
concept used in this appendix is more likely to furnish useful results 


52 


NACA TN 2608 


for stiffness distributions than for the root stiffnesses, and, because 
of the type of structure assumed, the concept is not applicable to very 
thin wings. 


Assumed Applied Loads 

If the applied normal load is distributed in a manner proportional 
to the chord, that is, if 

l = Kc (Bl) 

the bending moment at any point on the span can be obtained by integrating 
the chord distribution as follows: 





ds* ds* 


where s* is the dimensionless distance along the reference axis 
measured from the effective root. Similarly, the total normal load on 
one wing is given by 



c ds* 


If the wing is linearly tapered, so that 


c = c r [l - (1 - X) s*| 
where the taper ratio X is defined by 


(B2) 


X = 


c t 




then the ratio of the bending moment at any point of the span to the 
product of the total normal load and the wing semispan less one-half of 
the fuselage width can be expressed as follows: 


P 


M 

b’/2 
cos A • 


f 6 (s*,X) 


(B3) 
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where the function fg of s* and X is defined by 

f 6 5 j(rfr - rrx ^ - s *> 2 <>*> 

and shown in figure 19. 

Similarly, if the moment arm of the normal load applied to the wing 
at any station is also proportional to the chord, the constant of propor- 
tionality being e^, the distributed torque at any station is then e^cZ, 

and the accumulated torque is 


T = e l K ls& jH (B5) 

% 

which may, in turn, be expressed as 


i ^ ? =f 7 (s*,X) - • (B 6) 


where the function fj of s* and X is defined by 


f -Jt 

7 " 3 


1 + X + (2 4 - A,) ( 1 - A,) 


(1 + X) 2 (1 + x) 2 

and the average chord c is defined by 


s* + 


(1 + X) 2 


(1 - s*) (B7) 


c = Cr + C t 
2 

The function fj is also shown in figure 20. 


(B8) 


The total nonnal load on one wing, P, can be estimated from the 
design gross weight and the design load factor of the airplane in the 
following manner: 


p = k L w - nW w) 


(B9) 
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If the fraction of the wing lift to the total lift carried by the airplane 
(including that of the fuselage and tail) is designated by t^, so that 


’ll = 


Lw 

^Total 


and the fraction of the wing weight (including the amount of fuel, 
external stores, and so on used in the critical design condition) to the 
total design gross weight is designated by t \ 2 > 


n - Ww 

n2 “T 


then equation (B9) may be written as 


P = i r| 3 nW 


where 


(BIO) 


^3 = ’ll - *1 2 

With the value of P given by equation (BIO), equations (B3) and (b 6) 
serve to express the local bending and torsional moments in terms of 
known design parameters; 


Effective Skin Thickness Required to 
Resist Applied Loads 

The wing structure has to resist both the applied bending moments 
and the applied torques) in other words, the load-carrying members must 
resist combined axial and shear stresses. A relation commonly used in 
the design of wing structures loaded by compressive and shear stresses 
due to bending and torsion moments is 
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where f^ is the applied bending stress, the applied shear stress, 

F-g the allowable (compressive) bending stress, and Frp the allowable 
shear stress. However, a somewhat conservative relation. 



(Bll) 


is much more convenient for the present purpose and, consequently, is 
used as the basis of the following development. If the margin of safety 
is not zero, equation (Bll) can be rewritten as 



(B12) 


where is an effectiveness factor which can be expressed in terms of 

the margin of safety (M.S.) as 


1 

% " 1 + M.S. 


(B13) ‘ 


The applied bending stress is 


Mz 

I 


(B14) 


where z is the maximum ordinate on the compression side measured from 
and normal to the chordwise principal axis. Similarly, the applied shear 
stress is 


f t “ 2 At ( B15 ^ 


where A is the cross-sectional area of the (assumed) single torsion cell, 
and t is the skin thickness on the compression side. Substitution of 
equations (Bl4) and (B15) into equation (B12) yields 


Mz L + T I 

rF B l M f t 2Atz 


= v 


(Bl6) 
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In order to relate the bending and torsion stiffness of the wing to 
the skin thickness t or to an equivalent thickness t e which includes 
the material in the stringers and spar flanges, the bending stress is 
assumed to be carried by a box covered with sheet of an effective thick- 
ness t , the webs are assumed to carry no bending stress, and the 

torsion stress is assumed to be resisted by an equivalent single cell 
the two webs of which contain all the material of the actual webs. The 
torsion and .bending moments of inertia may then be written as 


and 




ch 2 t 

e 


1 = V‘e I 




11 


2 ^ 2„ 2 /h 

n 12 n l3 (2 


(B17) 


■ T) 


15 



(B18) 


where the effectiveness factors to are defined in table 1. 

In the factor t^, the effective perimeter p. of the torsion cell is 

the sum of the lengths of skin around the perimeter, each weighted by the 
ratio of the thickness of the critically stressed element to the thickness 
of the given length of skin. 

When the value of I given by equation (Bl8) is substituted into 
equation (Bl6) , equation (Bl6) may be written as 


^e 



n 15 c 



2Azt] 6 


(B19) 


By making use of equations (B3), (B5), and (BIO), as well as the effective- 
ness factors T)^ to tuq defined in table 1, equation (B19) can be 
written as y 
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t e = 


nW b'/2 


T l4 T li5 f B c 2 cos A h / c 


^i^i _i_ 
a a f 8 


The factor f8 is defined in terms of the factors fg and f^ given 
by equations (B4) and (B7) as 


f = - — 

8 2 tj 


(B21) 


The function fo is shown in figure 21. 


Bending and Torsional Stiffnesses 

Substitution of the value of t e given by equation (B20) into 
equation (Bl8) yields an expression for the bending moment of inertia I 
or for the bending stiffness El at any point along the span. The value 
of this stiffness at any point on the span may be divided by the stiff- 
ness (El) r at the wing root. This ratio can then be expressed as 


f + Vi 

_h_ £6_ _j3_ , 8 A A 

f 6 r f 8 r , . ¥l 

■ 8 r Aa 


... h/c /* x ^19 e l 

(h/c)„ f 9\ s a a , 


(B22) 


where 


Tl l9 e l 


f c f 7, 8 V 
9 " C r f ^r f Q + n 19 e r 

r Ax 


(B23)' 


The function f^ is plotted in figure 22. The value (El) r may be 
obtained from equations (Bl8) and (B20) as 
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<“>r -30 f 


where 


= (1 + ^) 3 %r( fQ + ¥l] 

r 64 f 8 r \ 8r a A j 


(B24). 


(B25) 


The function F r is shown in figure 3 as a function of X, 
as the parameter. 


with 

a A 


Similarly the torsional stiffness GJ may be obtained by substituting 
the value of tg given by equation (B20) into equation (B17). However, 
from equations (B17) and (Bl8) the ratio of the torsional stiffness to 
the bending stiffness may be obtained in the form 


GJ = 8 fi (B2 6) 

El E T] 8 T l 9 T ll5 


This equation shows that the ratio GJ/EI is constant along the span 
within the framework of the constant-stress concept • Equation (B26) may, 
therefore, he interpreted as an expression for the value GJ/EI at the 
wing root, that is, for the value (GJ) r /(El) r . The torsional stiffness 
at any other point on the span can then he . obtained from equation (B22), 
since 


GJ KE 

tgj!7 wt; 


(B27) 


because GJ/EI is constant over the span. 


The stiffness ratios GJ/(GJ) r and El/(El) r can be obtained 
directly from figure 22 when the thickness ratio h/c of the wing is 
constant along the span) if the thickness ratio is not constant the 
factor fo obtained from figure 22 must be multiplied by the ratio 
h/ c 

^ h/c ) " ' a "k an J station to obtain the stiffness ratio at that station. 

As may be seen from figure 22, the function f^ does not vary much with 

Wi. 


the parameter 


-j -this parameter represents the additional amount of 


A 


NACA TN 2608 


59 


skin thickness required to carry the torque (see equations (B19) 

and (B20)), and this additional thickness is small for most conventional 


wing structures. Consequently, an average value of 


Mi = 0.03 
a a 


was 


used to obtain the stiffness distributions used in the aeroelastic 
calculations on which the charts of this paper are based. 


Equation (B22) shows that, once a value has been assumed for the 


term 


^I9 e l 


a a ' 


the stiffness ratios El/(El) r and GJ/(GJ) r are inde- 


pendent of the effectiveness factors used in this analysis. Therefore, 
specific values of these parameters need not be known in order to estimate 
the stiffness distributions, but one of the assumptions on which equa- 
tion (B22) is based is that whatever values the effectiveness factors 
have are nearly constant along the span. In order to estimate the value 
of (El) r , however, these factors must be known, since they enter 
directly into equation (B24) . The estimate of (El) r obtained in this 
manner is, therefore, subject to all the limitations imposed by the 
approximations of the const ant -stress concept. Hence, some judgement 
must be, exercised in using this estimate, and, if possible, it should be 
modified in the light of . experience . 


Structural Weight Associated with the Stiffness Distribution 

The increase in structural weight associated with a given increase 
in stiffness can be estimated on the basis of assumptions similar to 
those made in relating the stiffness and the strength. For the purpose 
of this analysis the various components of the wing structure, exclusive 
of the carry-through structure within the fuselage, are classified in 
two groups, one which contains the elements that take the bending and 
torsional loads due to the assumed loading and one which contains all 
other components. In the first group are: 

(1) The amount of top and bottom skin that is used in the estimation 
of the thicknesses required to withstand the bending and torsional loads, 
including stringers and spar flanges included in the equivalent skin 

(2) Webs, including any web stiffeners 
In the second group are: 

(1) Skin, stiffeners, false spars, and so on, which are not con- 
sidered in the estimation of the equivalent thicknesses 

(2) Ribs, bulkheads, and posts designed to raise the buckling 
strength of the cover sheets 



60 


NACA TN 2608 


(3) Control surfaces and their supports, attachments, and actuating 
mechanisms 

( h) The supports of internal stores 

This analysis is concerned only with the first group and, more 
specifically, with the relative increase in the weight of this group 
occasioned by an increase in stiffness of the main structure. Means of 
estimating the actual magnitude of the weights involved and of est imat ing 
the weights of the items in the second group as well are given in refer- 
ences 7 and 8. 

The weight per unit length of the structural elements of the first 
group can be written as 


w s = 2l l8 T l2l’'8 ct e 


(B28) 


where 7 S is the density of the material of the primary structure (or 
an equivalent density in the case of sandwich construction), t) 21 is 
the ratio of an equivalent perimeter p to the actual perimeter of the 
cell, and p is the sum of all the lengths which constitute the perimeter, 
each multiplied by the ratio of its equivalent thickness to the equivalent 
thickness t e of the upper cover sheet. 


In view of the assumption made concerning the combination of bending 
and torsional stresses, the thickness t e required in equation (B28) can 
be obtained from equation (Bl8) as 


so that 


kl 

n 15 ch 2 




- 87, 


^ 8^21 _I_ 

^15 h2 


or 

_ f/fr 

Ws r ~ (h/h r ) 2 


(B29) 
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Consequently, the total weight (for both wings) of the structural elements 
of the first group can be estimated from the relation 


W s /2 


V 


s 


V/2 
r cos A 



as* 


(h/h r ) 2 


(B30) 


Equation (B30) serves to estimate relative changes in the weight of 
the first group of structural elements. For instance } with a given 
distribution of I and h^ that weight is directly proportional to I 
and inversely proportional to h^ 2 . Similarly, given two different P 

distributions of I r and h r with the same values at the root, the 
ratio of the weights is equal to the ratio of the two values obtained by 
using the respective distributions of I and h in the integral of 
equation (B30) . 

Although the actual value of W s is not relevant to this discussion, 
it may be estimated by substituting the previously calculated stiffness 
distributions into equation (B30), and the result is given here as a 
matter of general interest: 


where 


W = 2 J^W^l _nW_ b'/2 

s \\ 5\8 f b h A cos a w 


(B31) 


■ w 


32 

(1 + ? 0 2 

4 



-i^-ds* 

(h/h r ) 2 


(B32) 


According to equation (B31), the structural weight is directly 
proportional to the design gross weight, load factor, swept-span aspect 
ratio, span, and density of the material of the primary structure and 
inversely proportional to the allowable stress and the wing thickness 
ratio. The dependence of the weight on the taper ratio (all other param- 
eters, notably the aspect ratio and span, are the same) is illustrated in 
figure 14 by a plot of the function F w against taper ratio for several 


values of the parameter 


^19 e l 


and of the ratio of the wing thickness 


ratios at the tip and at the root. 
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TABLE 2 
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TABLE 3 .- PARAMETERS OF WING USED IN ILLUSTRATIVE EXAMPLE 


Geometrical parameters 

A 

k 

A, deg 

• 37.5 

S, sq in 

37,498 

b, in 

• 387.4 

v/ 2 , in 

. 20.0 

b'/ 2 , in 

• 173.7 

c r > ln ? 

. 102.8 

c t> in ' 

. 54.2 


• 0.527 

c MAC> ln 

. 100.4 


Structural parameters 


e 

. . o.kkk 

(GJ) r , lb/ sq in 

8.94 x io9 

(El) r , lb/ sq in 

9.56 x io9 

Gj/(GJ) r 

% (c/c r )^ 

El/(El) r > 

* (c/c r )^ 
. . 0 

As^, in 

As^, in 

. . '-5.6 

As, in 

. .. -3.0 

s t (.= ^-^-7 + ^ 

1 cos A J 

. - 215.9 


Aerodynamic parameters 


Subsonic 

Supersonic 


(M < 0.65) 

(M = 1.5) 

a 

0.25 

0.425 

e i 

0.194 

0.019 

lj ae 

2.78 

4.92 

K . . . 

0.78 

1.00 


" Aeroelastic parameters 


Subsonic 

Supersonic 


(M < 0.65) 

(M = 1.5) 

k 

7.76 

79.0 

K i 

2.82 

2.82 

K 2 

0.474 

0.474 

Q*D 

"1.053 

-0.0774 

q D , lb/ sq ft 

-6400 

. -2700 

f i 

Fig. 5(c) 

Fig. 5(c) . 

^1 

Fig. 5(c) 

Fig. 5(c) 

Fl 

-0.27 

-0.02 

V 1 • 

0.655 

0.662 

“l 

1.303 

1.333 

c w 

1 - 0.345=5- 
q D 

1 - 0 . 338 ^ 

CL w 0 

1 - .a. 

1 - -a. 


q D 

q D 


1 - 0.146=5. 

1 - 0.118=3- 

s 

q D 

q D 

*0 

. 1 - 0.345^ 

1 - 0.338jS. 

q D 
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Figure 2.- Charts for a preliminary survey of aeroelastic phenomena. 
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(c) Wings of taper ratio 0.5. 
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(e) Wings with moment arm ej_ = 0 (or | k| > 25). 


Figure 2. - Concluded. 
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(a) Effect of taper ratio on q*p for unswept wings. (q*^ = Kp in this 

case.) 

Figure U.- Effect of taper ratio on the dynamic pressure at divergence. 
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(b) Effect of taper ratio on q-p for e^ - 0. 



in this case 


) 


Figure h.~ Concluded. 
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(a) Stiffnesses given by constant-stress criterion for 


(h/c) t 

(h/c) r 


Figure 5>.- The angle-of-attack distribution functions f]_, Af]_, and 

F]_ for constant geometric angles of attack. t 



NACA TN 2608 


i m 


Function ft 


0 4 

Parameter k 













V 

7 " 

s 

s / 






\F unction A {r 


(ptck 

(h/c) r 1 

t.O ( excess strength 


Dimensionless distance along the -span, s 


(b) Stiffnesses related to those given by constant-stress criterion for 

wings with taper ratio 0.$. 

Figure' Continued. 
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(c) Stiffnesses proportional to 
Figure Concluded. 
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(a) Stiffnesses proportional (b) Stiffnesses given by 

, t . stress criterion 

to c 4 . 

Figure 6.- The angle-of-attack distributions for constant ge 
of attack for wings of taper ratio 0. 
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(c) Stiffnesses proportional to c^. 


Figure 7 .- Concluded. 
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(a) Stiffnesses proportional ' (b) Stiffnesses given by constant- 

t 0 C U # stress criterion. 

Figure 8.- The angle-of-attack distributions for linearly varying angles 
of attack for wings of taper ratio 0. 




Asymptotes 
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(b) Stiffnesses related to those given by constant-stress ci 

wings with taper ratio 0.5. 


Figure 9 Continued. 


Asymptotes 
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(c) Stiffnesses proportional to 
Figure Concluded. 
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t> .5 


' C &0 -8 


. C 


r 0 -8 



f ■ 

Dynamic pressure ratio , 


(a) Stiffnesses proportional (b) Stiffnesses given by constant- 

, L stress criterion, 

to c u . 

Figure 10.- The lift- and moment-coefficient ratios for constant geometric 
angles of attack for wings of taper ratio 0. 
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(b) Stiffnesses related to those given by constant-stress criterion for 

wings with taper ratio 0.$. 


Figure 11.- Continued 


Asym+otes 
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A 

AO 

.5 



Parameter k 


(c) Stiffnesses proportional to 
Figure 11.- Concluded. 



Dynamic-pressure ratio , 


J±- 

(a) Stiffnesses proportional (b) Stiffnesses given by constant- 
bo c^. stress criterion. 


Figure 12.- The lift- and moment-coefficient ratios for linearly, varying 
geometric angles of attack for wings of taper ratio 0. 
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Figure 13.- Effect of aeroelastic action on some aerodynamic properties 
of the wing used in illustrative example . 
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Figure 16.- Comparison of angle-of-attack ratios calculated by the matrix 
method of appendix A with those calculated from equation (20) for 
constant geometric angles of attack at various dynamic -pres sure ratios 
X = 0.5, k = 8, stiffnesses proportional to c^. 
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Figure 19.- The bending-moment function f^ 



Figure 20.- The twisting -moment function f. 
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